Simulation of anyons with tensor network algorithms 
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Interacting systems of anyons pose a unique ciiallenge to condensed matter simulations due to 
their non-trivial exchange statistics. These systems are of great interest as they have the poten- 
tial for robust universal quantum computation, but numerical tools for studying them are as yet 
limited. We show how existing tensor network algorithms may be adapted for use with systems of 
anyons, and demonstrate this process for the 1-D Multi-scale Entanglement Renormalisation Ansatz 
(MERA). We apply the MERA to infinite chains of interacting Fibonacci anyons, computing their 
scaling dimensions and local scaling operators. The scaling dimensions obtained are seen to be 
in agreement with conformal field theory. The techniques developed are applicable to any tensor 
network algorithm, and the ability to adapt these ansatze for use on anyonic systems opens the door 
for numerical simulation of large systems of free and interacting anyons in one and two dimensions. 

PACS numbers: 05.30.Pr, 73.43.Lp, 02.70.-c, 03.65.Vf 
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I. INTRODUCTION 

The study of anyons offers one of the most exciting 
challenges in contemporary physics. Anyons are ex- 
otic quasiparticles with non-trivial exchange statistics, 
which makes them difficult to simulate. However, they 
are of great interest as some species offer the prospect 
of a highly fault-tolerant form of universal quantum 
computational^ and it has been suggested^ that the sim- 
plest such species may appear in the fractional quantum 
Hall state with filling fraction v = 12/5. Despite the 
current strong interest in the development of practical 
quantum computing, our ability to study the collective 
behaviour of systems of anyons remains limited. 

The study of interacting systems o f anyons usin g 
numerical techniques was pioneered by iFeiguin et al\ A 
using exact diagonalisation for 1-D systems of up to 
37 anyons, and the Density Matrix Renormalisation 
Group algorith m (DMRGyS for long er chains. Also re- 
lated is w ork by ISierra and Nishinofe later extended by 
iTatsuakl — which applies a variant of DMRG to spin 
chain models having SU{2)k symmetry. Some of these 
models are now known to correspond to SU{2)k anyon 
chains,^ and using this mapping these systems may also 
be studied using the Bethe ansatz^ and quantum Monte 
Carlo.ift. 

However, all of these methods have their limitations. 
Exact diagonalisation has a computational cost which is 
exponential in the number of sites, strongly limiting the 
size of the systems which may be studied. DMRG is 
capable of studying larger system sizes, but is typically 
limited to 1-D or quasi-l-D systems (e.g. ladders). Map- 
ping to a spin chain is useful in one dimension but is 
substantially less practical in two. There are therefore 
good reasons to desire a formalism which will allow the 
application of other tensor network algorithms to sys- 
tems of anyons. Many of these tensor networks, such as 
Projected Entangled Pair States (PEPS)fii"i^ and the 



2-D versions of Tree Tensor Networks (TTN)^'^ and of 
the Multi-scale Entanglement Renormalisation Ansatz 
(MERA)i2i"— have been designed specifically to accu- 
rately describe two-dimensional systems. 

In one dimension, many previously studied systems 
of interacting anyons display extended critical phases;^ 
which are characterised by correlators exhibiting poly- 
nomial decayj^ Whereas DMRG favours accurate rep- 
resentation of short range correlators at the expense of 
long-range accuracy, the 1-D MERA^^'^^ is ideally suited 
to this situation as its hierarchical structure naturally 
encodes the renormalisation group flow at the level of 
operators and wavefunctions,^^^^"* and hence accurately 
reproduces correlators across a wide range of length 
scalesj ^^i^^'^^" — The development of a general formal- 
ism for anyonic tensor networks is therefore also advan- 
tageous for the study of 1-D anyonic systems. 

This paper describes how any tensor network algo- 
rithm may be adapted to systems of anyons in one or two 
dimensions using structures which explicitly implement 
the quantum group symmetry of the anyon model. As a 
specific example we demonstrate the construction of the 
anyonic 1-D MERA, which we then apply to an infinite 
chain of interacting Fibonacci anyons at criticality. The 
approach which we present is completely general, and 
can be applied to any species of anyons and any tensor 
network ansatz. 



II. ANYONIC STATES 

Consider a lattice Cq of n sites populated by anyons. 
In contrast to bosonic and fermionic systems, for many 
anyon models the total Hilbert space Yco can not be 
divided into a tensor product of local Hilbert spaces. In- 
stead, a basis is defined by introducing a specific fusion 
tree (e.g. Fig.[T][i)). The fusion tree is always constructed 
on a linear ordering of anyons, and while the 1-D lattice 
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FIG. 1. (i) Example representation of a state [i/') in a fusion 
tree basis for a system of 6 anyons. Labels ai indicate charges 
associated with edges of the fusion tree graph, and labels Ui 
are degeneracies associated with vertices. The structure of 
the tree corresponds to a choice of basis, and does not affect 
the physical content of the theory, (ii) Braiding may be used 
to change the ordering of the leaves of a fusion tree basis, or 
to represent anyon exchange, (iii) _F-moves convert between 
the bases associated with different fusion trees. 

naturally exhibits such an ordering, for 2-D lattices some 
linear ordering must be imposed. Each line is then la- 
belled with a charge index ai such that the labels are 
consistent with the fusion rules of the anyon model, 

axb^Y.^abC- (1) 

c 

For anyon types where some entries of the multiplicity 
tensor A^^^^ take values greater than 1, a label Ui is also 
affixed to the vertex which represents the fusion process 
to distinguish between the different copies of charge c. 
The edges of the graph which are connected to a vertex 
only at their lower end are termed "leaves" of the fusion 
tree, and we will associate these leaves with the charge 
labels ai...a„. Different orderings of the leaves on a 
fusion tree may be interconverted by means of braiding 
(Fig. HJii)), and different fusion trees, corresponding to 
different bases of states, may be interconverted by means 
of F moves (Fig. [ljiii))i22,i^ In some situations it may 
also be useful to associate a further index bi with each of 
the leaves of the fusion tree. For example, if the leaves 
are equated with the sites of a physical lattice, then this 
additional index may be used to enumerate additional 
non-anyonic degrees of freedom associated with that lat- 
tice. For simplicity we will usually leave these extra in- 
dices hi . . . bn implicit, as we have done in Fig. [TJ as they 
do not directly participate in anyonic manipulations such 
as F moves and braiding. 

Let the total number of charge labels on the fusion 
tree be given by m, where m > n. For abelian anyons 
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FIG. 2. The leaves of this fusion tree carry the charge labels ai 
to ag. An edge which is not a leaf, labelled with charge aio, is 
indicated by the large grey arrow. The portion of the fusion 
tree extending from edge aio out to the leaves is indicated 
by the grey ellipse. If a degeneracy index /iio is associated 
with charge aio, then for a given value of aio, index /xio will 
enumerate all compatible labellings of the highlighted portion 
of the fusion tree. 

the fusion rules uniquely constrain all a.i for i > n, and 
provided there are no constraints on the total charge, the 
total Hilbert space reduces to a product of local Hilbert 
spaces V, such that V^^ = V®". For nonabelian anyons, 
additional degrees of freedom arise because some fusion 
rules admit multiple outcomes, permitting certain ai {i > 
n) to take on multiple values while remaining consistent 
with the fusion rules, and the resulting Hilbert space does 
not necessarily admit a tensor product structure. 

We will now associate a parameter Vi^ai with each 
charge on the fusion tree, which we will term the de- 
generacy. This parameter corresponds to the number of 
possible fusion processes by which charge ai may be ob- 
tained at location i. Where charge a^ arises from the 
fusion of charges and aj, then i^k.ak will satisfy 

•^k^a, = ''^,a,''J,a,K'ar (2) 

ai ,aj 

For systems where the only degrees of freedom are any- 
onic, degeneracies on the physical lattice Co (i.e. Vi^ca, 
I < i < n) will take values of or 1 depending on whether 
a charge ai is permitted on lattice site i. Higher values 
of Ui^at may be used on the physical lattice if there is 
also a need to represent additional non-anyonic degrees 
of freedom, enumerated by indices bi . . .bn- 

Up to this point we have parameterised our Hilbert 
space in terms of explicit labellings of the fusion tree. 
We now adopt a different approach: Consider an edge 
i of the fusion tree which is not a "leaf" . As well as 
labelling this edge with a charge ai we may introduce a 
second index fj,i, running from 1 to i^i,ai- Each pair of 
values {aj,/ij} may be associated with a unique charge 
labelling for the portion of the fusion tree from edge i 
out to the leaves, with these labellings being compatible 
with the fusion rules in the presence of a charge of ai on 
site i (for an illustration of this, see Fig. [5]). Provided we 
know the structure of the fusion tree above i and have a 
systematic means of associating labellings of that portion 
of the tree with values of ^i , then in lieu of stating the 
values of all aj for edges j involved in that portion of the 
tree, we may simply specify the value of the degeneracy 
index ^i. In this way we may specify an entire state in 



3 



T T T T 



f =1, 



T T T T 








FIG. 3. Example enumeration of states according to am and 
^rn for a fusion tree describing four Fibonacci anyons. The 
Fibonacci anyon model has one non-vacuum charge label (r) 
and one non-trivial fusion rule, r x r — >■ 1 + r. Because the 
charges 1 and r are both self-dual, no arrows are required 
on diagrammatic representations of Fibonacci anyon fusion 
trees. 



the form 



1^) 



Mm 



1 1 f^m ) 



(3) 



where am is the total charge obtained on fusing all the 
anyons. The index /i^, which is the degeneracy index 
associated with the total charge of the fusion tree, may 
be understood as systematically enumerating all possi- 
ble labellings of the entire fusion tree including charge 
labels, vertex labels, and any labels associated with ad- 
ditional non-anyonic degrees of freedom. For an example, 
see Fig. [3l Note that for a given edge i, the value of the 
degeneracy Vi^ai may vary with the charge Oi and conse- 
quently the range of the degeneracy index fim in Eq. ([3]) 
is dependent on the value of the charge am- 

The notation of Eq. ([3]) should be contrasted with that 
of Fig. [iji). In the latter, the number of indices on c 
depends upon the number of charge labels on the fu- 
sion tree, whereas in the former, the tensor describing 
the state is always indexed by just one pair of labels — 
charge and degeneracy — which will prove advantageous 
in constructing a tensor network formalism for systems 
of anyons. 

We now choose to restrict our attention to systems hav- 
ing the identity charge. We may do this without loss of 
generality as a state on n lattice sites with a total charge 
a„i may always be equivalently represented by a state on 
n + 1 lattice sites whose total charge is the identity, with 
a charge on lattice site n + 1. This additional charge 
annihilates the total charge Um of sites 1 ... n to give the 
vacuum. The expression for 1-0) then becomes 



|l,Mm') 



(4) 



where fim' ranges from 1 to the dimension of the Hilbert 
space of the system of n sites with total charge am- Con- 
sequently we may represent the state ji/i) of a system of 



anyons by means of the vector ci^^, . For simplicity of 
notation, we will take greek indices from the beginning 
of the alphabet to correspond to pairs of indices {at, ^i} 
consisting of a charge index and the associated degener- 
acy index. The vector ci^^, will therefore be denoted 
simply c", with the understanding that in this case the 
charge component am' of multi-index a takes only the 
value 1. (Multi-index a is raised as we will shortly in- 
troduce a diagrammatic formalism in which vector c is 
represented by an object with a single upward-going leg. 
In this formalism, upward- and downward-going legs may 
be associated with upper and lower multi-indices respec- 
tively.) 



III. ANYONIC OPERATORS 

We will divide our consideration of anyonic opera- 
tors into two parts. First we shall consider operators 
which map a state on some Hilbert space H into an- 
other state on the same Hilbert space. When applied to 
a state represented by c", such an operator leaves the 
degeneracies of the charges in multi-index a unchanged. 
We will therefore call these degeneracy-preserving any- 
onic operators. Then we will consider those operators 
which map a state on some Hilbert space T-L into a state 
on some other Hilbert space H' - These operators may 
represent processes which modify the environment, for 
example by adding or removing lattice sites, and also 
play an important part in anyonic tensor networks, for 
instance taking the role of isometrics in the TTN and 
MERA. As these operators can change the degenera- 
cies of charges in a multi-index a, we will call them 
degeneracy-changing anyonic operators. More generally, 
the degeneracy-preserving anyonic operators may be con- 
sidered a subclass of the degeneracy-changing anyonic 
operators for which T-L = %' . 



A. Degeneracy preserving anyonic operators 

We begin with those operators which map states on 
some Hilbert space H into other states on the same 
Hilbert space H. Examples of these operators include 
Hamiltonians, reduced density matrices, and unitary 
transformations such as the disentanglers of the MERA. 

First, we introduce splitting trees. The space of split- 
ting trees is dual to the space of fusion trees. While 
the space of fusion trees consists of labelled directed 
graphs whose number of branches increases monotoni- 
cally when read from bottom to top, the space of splitting 
trees consists of labelled directed graphs whose number 
of branches increases monotonically when read from top 
to bottom. An inner product is defined by connecting the 
leaves of fusion and splitting trees which have equivalent 
linear orderings of the leaves (braiding first if necessary) , 
then eliminating all loops as per Fig. |4l[i), with F moves 
performed as required. 
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FIG. 4. (i) Loops are eliminated by replacing them with 
an equivalent numerical factor determined by the normalisa- 
tion convention. The factor given here corresponds to the 
diagrammatic isotopy convention employed in Ref [2^ . (ii) 
Definition of a simple two-site anyonic operator, (iii) Appli- 
cation of an operator to a state is performed by connecting 
the diagrams' free legs. By performing F moves and elimi- 
nating loops (and in more complex examples, also braiding) 
it is possible to obtain an expression for the resulting state in 
the original basis. 



FIG. 5. An operator acting on four Fibonacci anyons. The 
values of the coefficients Matcde rnay be specified as a block- 
diagonal matrix Mcf , for example as in Table |ll 

Similarly, values of /3 correspond to pairs {aj,fJ-j} where 
flj has degeneracy . For degeneracy-preserving any- 
onic operators the charge indices and aj necessarily 
take on the same range of values, and z/q. = Va when 
fli aj. The values of may equivalently be calcu- 
lated from either the fusion tree making up the top half 
or the splitting tree making up the bottom half of the 
operator diagram. 

A well-defined anyonic operator M must respect the 
(quantum) symmetry group of the anyon model, and con- 
sequently all entries in M^^ for which ^ aj will be 
zero. However, in contrast with c" we do not require 
that Qi = Qj = 1. When M is a degeneracy-preserving 
operator, matrix M^^ is therefore a square matrix of side 
length 

^M^^i'ao (5) 

which may be organised to exhibit a structure which is 
block diagonal in the charge indices and aj, and for 
which the blocks are also square. As an example consider 
Fig. [HI which shows an operator acting on four Fibonacci 
anyons. An example matrix M^^ for an operator of this 
form is given in Table |T1 from which the entries of Matcde 
can be reconstructed, e.g. Mririr = 3. 



Anyonic operators may always be written as a sum 
over fusion and splitting trees, such as the two-site oper- 
ator M shown in Fig.|4l[ii), and for degeneracy-preserving 
anyonic operators it is always possible to choose the split- 
ting tree to be the adjoint of the fusion tree. To apply an 
operator to a state the two corresponding diagrammatic 
representations are connected as shown in Fig.|4ljiii), and 
closed loops may be eliminated as shown in Fig.Ufi). Se- 
quences of F moves, braiding, and loop eliminations may 
be performed until the diagram has been reduced once 
more to a fusion tree without loops on a lattice of n sites. 

Much as the state of an anyonic system may be repre- 
sented by a vector c", anyonic operators may be repre- 
sented by a matrix AI^ ^ . Each value of a corresponds to 
a pair {ai,/ii} where at is a possible charge of the cen- 
tral edge of the operator diagram (e.g. as in Fig. IDJii)), 
and fj,i is a value of the degeneracy index associated with 
charge a^. We will denote the degeneracy of by Vai- 



B. Degeneracy changing anyonic operators 

We now introduce the second class of anyonic oper- 
ators, which map states in some Hilbert space T-L into 
some other Hilbert space H'. These operators may re- 
duce or increase the degeneracy of any charge present in 
the spaces on which they act, and may even project out 
entire charge sectors by setting their degeneracy to zero. 
When these operators are written in the conventional no- 
tation of Fig. 21 the fusion and splitting trees will not 
be identical. Further, we may choose to allow combina- 
tions of degeneracies which do not naturally admit com- 
plete decomposition into individual anyons. For exam- 
ple, a degeneracy-changing operator may map a state on 
five Fibonacci anyons (having total degeneracies vi — 3, 
Vr — 5) into a state having degeneracies vi — 2, Vt — 2. 
As these degeneracies do not admit decomposition into 
an integer number of nondegenerate anyons, it is neces- 
sary to associate an index ui with the single open leg of 
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TABLE I. Matrix representation MJ' for an example opera- 
tor of the form shown in Fig. [S] Multi-index a corresponds 
to index pair {ai,/^^} and multi-index /? corresponds to pair 
{a-i, /Xj}. Subject to an appropriate ordering convention for 
\ii and these indices may be related to the fusion tree la- 
bels o, b, c, d, e of Fig.Oas shown. Note that as c is the charge 
on the central leg of Fig. [5] all nonzero entries of MJ' satisfy 



fli = a, = c. 



the fusion tree. This index behaves identically to the ver- 
tex indices Ui of Fig. [1] serving to enumerate the different 
copies of each individual charge, and as with the vertex 
indices of Fig. [1] it is absorbed into the degeneracy index 

As a further example, a state having degeneracies v\ = 
4, i/T- = 4 could be associated with a fusion tree having 
either one leg, or two legs each with degeneracies vi = 0, 
Vj- = 2. Again, indices Ui would have to be associated 
with each open leg. 

Matrix representations of degeneracy-changing any- 
onic operators may also be constructed, and when they 
are written in block diagonal form, the matrices and 
their blocks may be rectangular rather than square. 
Degeneracy-changing anyonic operators therefore repre- 
sent a generalisation of the degeneracy-preserving any- 
onic operators discussed in Sec. IIII Al It is worth not- 
ing that the presence of indices Ui on the open legs of 
the fusion or splitting trees of an operator do not auto- 
matically imply that it is a degeneracy-changing anyonic 
operator: The defining characteristic of a degeneracy- 
preserving anyonic operator is that it maps a state in a 
Hilbert space T-L into a state in the same Hilbert space 
and consequently both the matrix as a whole and all 
of its blocks are square. Thus a degeneracy-preserving 
anyonic operator may act on states having additional in- 
dices Ui on their open legs, and the resulting state may 
be expressed in the form of the same fusion tree, with 
the same additional indices on the open legs. 

Operators which change degeneracies may represent 
physical processes which change the accessible Hilbert 
space of a system. As we will see in Sec. IV Al they may 
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FIG. 6. Diagrammatic representation of operators P*^^^ (|6]) 
and P^^^ (O. The charges on all leaves are non-degenerate. 



also be used in tensor network algorithms as part of an 
efficient representation of particular states or subspaces 
of a Hilbert space, for example the ground state or the 
low energy sector of a local Hamiltonian. 

This distinction between degeneracy-changing and 
degeneracy-preserving anyonic operators is clearly seen 
with a simple example. Let be a state on six Fi- 
bonacci anyons. This state can be parameterised by a 
vector c", which has five components. We now define 
two projection operators, P'^-' and P'-^' (Fig. E]), each 
of which acts on the fusion space of anyons ri and T2. 
Operator P^^^ is degeneracy-preserving, and projects c" 
into the subspace in which anyons ri and T2 fuse to the 
identity. Its matrix representation is 



,(1) _ 




(6) 



where the first value of each multi-index corresponds to 
a charge of 1, and the second to a charge of t. Operator 
p(2) performs the same projection, but is degeneracy- 
changing. Its matrix representation is written 



P 



(2) P 



(10). 



(7) 



Both operators perform equivalent projections, in the 
sense that 



(^|p(l)tp(l)|^) = (V'|P(2^tp(2)|^^)^ 



(8) 



When P*-^^ acts on it leaves the Hilbert space un- 
changed, and hence the vector c'" describing state — 
P'-^^|'0) is once again a five-component vector, although 
in an appropriate basis some components will now nec- 
essarily be zero. In contrast P*^^-' explicitly reduces the 
dimension of the Hilbert space, and the vector c"" de- 
scribing state li/)") = P^^'ji/i) is of length two, describ- 
ing a fusion tree on only four Fibonacci anyons (as both 
Ti and T2 have been eliminated). One consequence of 
this distinction is that while (P^^^)^ — P^^\ the value of 
(p(2))2 undefined. 



IV. ANYONIC TENSOR NETWORKS 
A. Diagrammatic notation 

The diagrammatic notation conventionally employed 
in the study of anyonic systems, and used here in Figs.[T] 
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FIG. 7. (i) Diagrammatic representation of a state and 
two-site operator M expressed in terms of degeneracy indices, 
(ii) Application of M to state \'ip). Grey shapes represent ten- 
sors with charge and degeneracy multi-indices, with each leg 
of the shape corresponding to one charge and degeneracy in- 
dex pair. These diagrams represent the same state, operator, 
and process as Fig. [TJi) and Fig. [4ljii)-(iii). 



and m is well suited to the complete description of any- 
onic systems, as it provides a physically meaningful de- 
piction of the entire Hilbert space. However, the number 
of parameters required for such a description grows expo- 
nentially in the system size, and because it is necessary to 
explicitly assign every index to a specific charge or degen- 
eracy, specification of a tensor network rapidly becomes 
inconveniently verbose (for example see Fig. |4l[iii)). 

In the preceding Sections, we developed techniques 
whereby anyonic states and operators could be repre- 
sented as vectors and matrices, bearing only one or two 
multi-indices apiece. We now introduce the graphical no- 
tation which complements this description, and in which 
we will formulate anyonic tensor networks. Fig.[7ji) gives 
the graphical representations of a state associated 
with a vector c", and of an operator M associated with 
a matrix M^^ . The circle marked c corresponds to the 
vector c", and the circle marked M corresponds to the 
matrix M^^ . In general, grey circles correspond to ten- 
sors, and the number of legs on the circle corresponds 
to the number of multi-indices on the associated tensor. 
Each multi-index is also associated with a fusion or split- 
ting tree structure, which is specified graphically. For 
reasons to be discussed shortly, we will require that no 
tensor ever have more than three multi-indices. As the 



legs of the grey shapes are each associated with a multi- 
index, they carry both degeneracy and charge indices. 
Consequently it is not necessary to explicitly assign la- 
bels to the fusion/splitting trees, as these labellings are 
contained implicitly in the degeneracy index (for example 
see Table U where specifying the values of {ai,^i} and 
{aj , fij } is equivalent to fully labelling the fusion and 
sphtting trees of Fig. [5]). 

The fusion or splitting tree associated with a particu- 
lar multi-index may be manipulated in the usual way by 
means of braids and F moves, recalling that each compo- 
nent of the tensor is associated with a particular labelling 
of the fusion and splitting trees via the corresponding val- 
ues of the multi-indices. Manipulations performed upon 
a particular tree thus generate unitary matrices which act 
upon the multi-index that corresponds to the labellings 
of that particular tree. 

The application of an operator to a state is, unsur- 
prisingly, performed by connecting the appropriate dia- 
grams, as shown in Fig. [TJii). For operators of the type 
discussed in Sec. IIII Al the outcome is necessarily a new 
state in the same Hilbert space, which consequently can 
be described by a new state vector c'", as shown. How- 
ever, in general an operator M will not act on the entire 
Hilbert space of the system, and so will be described by 
a tensor constructed on the fusion space of some subset 
of lattice sites, and not on the system as a whole. Oper- 
ator M acting on state jV') in Fig. [7{ii) is an example of 
this. Because c" describes a six-site system but M^^ is 
constructed on the fusion space of two sites, the multi- 
indices of ^ span a significantly smaller Hilbert space 
than that of c" and we cannot simply write 

c"5 = c'^mJ (9) 

(using Einstein notation, where repeated multi-indices 
are assumed to be summed). Instead, we must under- 
stand how to expand the matrix representation of an op- 
erator on some number of sites x, to obtain its matrix 
representation as an operator on x' sites, where x' > x. 



B. Site expansion of anyonic operators 

The multiplicity tensor N^f^ describes the fusion of two 
charges without degeneracies. It is easily extended to 
incorporate degeneracies of the charges, and we will de- 
note this expanded multiplicity tensor A'^^^ where multi- 
indices a, /3, and 7 are associated with the pairs {a, /ia}, 
{6, /ifc}, and {c, /id respectively, and for given values of 
a, /3, and 7, u runs from 1 to A^^f,. The degeneracies asso- 
ciated with charges a, b, and c are denoted Va, i^b, and Vc 
respectively. As with fia, fJ-b, and fj,c, there is an implicit 
additional index on each degeneracy i^x representing the 
edge of the tree on which charge x resides. The values of 

and Vh may be chosen arbitrarily (for example, h'a\a=i 
may differ from Vb\b=i), but the degeneracies associated 
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Pair 


Assigned pentuplet 


{c,fic} 


{a, ^a,b, ^lb,u} 


1, 1 


1, 1, 1, 1, 1 


1, 2 


r, 1, r, 1, 1 


1, 3 


r, 1, r, 2, 1 


r, 1 


1, 1, r, 1, 1 


T, 2 


1, 1, r, 2, 1 


T, 3 


r, 1, 1, 1, 1 


T, 4 


r, 1, r, 1, 1 


r, 5 


r, 1, r, 2, 1 



TABLE II. Construction of N^^^ for a fusion vertex for Fi- 
bonacci anyons. In this example a may take charges 1 and 
T each with degeneracy 1, and b may take charges 1 and r 
with degeneracies 1 and 2 respectively. By Eq. ([2]), charge c 
may therefore take values 1 and r with degeneracies 3 and 5 
respectively. A correspondence between the values of multi- 
index 7 and of multi-indices a and /3 is established in some 
systematic manner, with each assignation satisfying c £ axb, 
and for Fibonacci anyons the index u is trivial as all multi- 
plicities N^fj are zero or one. An example assignation is shown 
in the table. The corresponding entries of are then set to 
1, with all other entries zero. For example, the fourth row 
indicates that iV/j'A', =1. 



with the values of c must satisfy 

,y, = Yl "-''bKb (10) 



in accordance with Eq. 1^. When this constraint is satis- 
fied, every quadruplet of indices {a, ^a, b, /Zb} correspond- 
ing to a unique pair of choices for a and /? may be asso- 
ciated with N^i^ distinct pairs of indices {c, ^c} for each 
c £ a X b. These pairs {c, /id are enumerated by the 
additional index u. This defines a 1:1 mapping between 
sets of values on {a, fia, b, fib, u} and pairs {c, fic}, and we 
set the corresponding entries in N^iSu 1' '^ith all other 
entries being zero. A simple example is given in Table HIl 



By virtue of their derivation from N^f^, the object 

7 



^apu ^^'^ conjugate N 



represent application of 
the anyonic fusion rules, and may be associated with 
vertices of the splitting and fusion trees. Under the iso- 
topy invariance convention there is an additional factor of 
[dc/ {dadb)]^ associated with the fusion of charges a and 
b into c, where dx is the quantum dimension of charge x, 
and similarly for splitting, but we will account for these 
factors separately. Thus constructed, the tensors N sat- 

When used as a representation of the fusion rules, the 
generalised multiplicity tensor N2p^ and its conjugate 
^'fafiu pgj-j^j^ increase or decrease the number of 

multi-indices on a tensor in a manner which is consistent 
with the fusion rules of the quantum symmetry group. 
This process is reversible provided the symmetry group 
is abelian or, for a nonabelian symmetry group, provided 



the total number of multi-indices on the tensor does not 
at any time exceed three. In constructing and manipu- 
lating a tensor network for a system of anyons, we will 
require only objects which respect the fusion rules of the 
anyon model. It is a defining property of such objects 
that when the number of multi-indices they possess is 
reduced to 1 by repeated application of N and N\ non- 
zero entries may be found only in the vacuum sector. We 
imposed this requirement for states in Sec. |lTl and it is 
equivalent to the restriction we imposed on anyonic op- 
erators in Sec. IIII Al In Ref. 30 an equivalent condition 
was observed for tensors remaining unchanged under the 
action of a Lie group, and these tensors were termed in- 
variant. When working with invariant tensors, we may 
separately evaluate the components of the tensors act- 
ing on the degeneracy spaces (e.g. the nonzero blocks 
of M^^), and the factors arising from loops and vertices 
of the associated spin network. This property greatly 
simplifies the contraction of pairs of tensors. 

In addition to increasing or decreasing the number of 
legs of a tensor, we may also use N to "raise" the matrix 
representation of an operator from the space of x sites to 
the space of {x + x') sites. This is shown in Fig. [51 and 
the matrix representation of the raised operator is given 
by 



Seu 



(11) 



where multi-index e describes the fusion space of all sites 
in (x + x') but not in x. Because the numeric factors 
associated with loops, vertices, and braiding (where ap- 
plicable) are handled separately, no factors of quantum 
dimensions appear in Eq. (ITT]) . 

To act an operator M on a state {ip) in the matrix 
representation, we therefore connect the diagrams for M 
and IVj), eliminate all loops, and then raise the matrix 
representation of the operator M using Eq. pT|) . repeat- 
edly if necessary, until the resulting matrix M'^ ^ may be 
applied directly to the state vector c" . Similarly it is pos- 
sible to combine the matrix representations of operators, 
by connecting their diagrams appropriately, eliminating 
loops, and performing any required raising so that both 
operators act on the same fusion space. Their matrix 
representations can then be combined to yield the ma- 
trix representation of the new operator: 



M 



(1x2) 



7^(2) ^ 



7 



(12) 



and the fusion/splitting tree associated with this new op- 
erator is obtained as shown in Fig. [S] 

Note that as yet, we have not described how two ob- 
jects may be combined if their multi-indices are both up 
or both down, and are connected by a curved line. To 
contract such objects together, it is necessary to under- 
stand how bends act on the central matrix of an opera- 
tor. Once this is understood, the bend can be absorbed 
into one of the central matrices, so that the connection 
is once again between an upper and a lower multi-index 
as in Eq. (fT^ . This process is described in Sec. IIV Cl 




FIG. 8. "Raising" of an operator M from sites x to sites 
X + x': (i) Operator M defined only on sites denoted x. (ii) 
Resolutions of the identity are inserted above and below M, 
being constructed from tensors N and N'' . The central por- 
tion of this diagram is identified as corresponding to the new 
matrix ^ which describes M on x + x'. (iii) Loop and ver- 
tex factors in the central region are evaluated separately and 
eliminated, (iv) The tensor network corresponding to the new 
central portion is contracted. The N and N"^ tensors outside 
the central region become vertices of the fusion and splitting 
trees associated with M^^ . Together the trees and the matrix 
M^'' constitute the raised version of M. 



C. Manipulation of anyonic operators 

As observed in Sec. IIVB[ when we describe a system 
entirely in terms of objects invariant under the action of 
the symmetry group, we may account separately for the 
numerical normalisation factors associated with the spin 
network. However, as well as affecting these numerical 
factors, transformations of the fusion or splitting tree of 
an anyonic operator will typically also generate unitary 
matrices which act on the matrix representation of the 
operator. These matrices respect the symmetry of the 
anyon model, and thus can be written as block-diagonal 
matrices where each block is a unitary matrix acting on 
a particular charge sector. In terms of the diagrammatic 
notation of Sec. HV^ F moves and braids therefore re- 
sult in the insertion of a unitary matrix, as shown in 
Fig. ini These matrices, whose entries are derived from 
the tensors {F^'"^)i^euv)(fu'v') and Rf' respectively, are 
raised if required, as described in Sec. IIVBI and then 
contracted with M^,^, the matrix representation of the 
operator. To compute the unitary matrices involved, it 
suffices to recognise that F moves and braids are unitary 
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FIG. 9. (i) F move, and (ii) braiding, performed on a section 
of fusion tree in the diagrammatic notation of Sec. IIV Al 



transformations in the space of labelled tree diagrams. 
Identifying the leg on which the unitary matrix is to be 
inserted, the relevant region of the space of labelled dia- 
grams is then enumerated by the multi-index which can 
be associated with this leg (compare Fig. [2]). 

Braiding is of particular importance when working in 
two dimensions, as an operator will necessarily be de- 
fined with respect to some arbitrary linear ordering of 
its legs, and when manipulating a tensor network it may 
be necessary to map between this original definition and 
other equivalent definitions, corresponding to different 
leg orderings. For example, let M be a four-site anyonic 
operator as shown in Fig. llOf i). which we wish to apply 
to a 2-D lattice. For the indicated linearisation of this 
lattice, application of M will require braiding as shown 
in Fig. llOf ii). By evaluating the unitary transformations 
corresponding to these braids and absorbing them into 
Mq,^, we may define a new operator M' which acts di- 
rectly on the linearised lattice without any intervening 
manipulations of the fusion/splitting trees. 

We will also frequently wish to deal with tensor legs 
which bend vertically through 180 degrees. If working 
with an anyon model that has non-trivial Frobenius- 
Schur indicators, then indicator flags must be applied 
to all bends. Like F moves and braiding, the reversal 
of a Frobenius-Schur indicator flag is a unitary transfor- 
mation, and once again this leads to the introduction of 
a unitary matrix which can be absorbed into a nearby 
existing tensor. However, we may wish to perform other 
operations on bends, such as absorbing them into fusion 
vertices or the central matrices of anyonic operators. We 
may also need to move a matrix ^ across a bend. We 
must therefore develop the description of bends in the 
new diagrammatic formalism. 

In Ref. 29 a prescription for absorbing bends into fu- 
sion vertices is given in terms of tensors {A'^'^)uv and 
{Bf )uv, derived from the F moves, and corresponding to 
clockwise and counter-clockwise bends respectively. The 
absorption of a clockwise or counterclockwise bend into 
a fusion vertex is reproduced in Fig. Illf i). and results in 
a vertex fusing upward- and downward-going legs. We 



(i) 
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(ii) 



M 
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13 2 4 



FIG. 10. (i) An operator M acting on sites on a 2-D lattice is 
defined with respect to some arbitrary linear ordering of these 
sites, (ii) When manipulating the tensor network, it may on 
occasion be computationally convenient for the lattice to be 
linearised according to some alternative linearisation scheme. 
In this example, the imposed linearisation scheme is indicated 
by the dotted line. To apply M to a different linearisation of 
the lattice may require braiding. The orientation of the braids 
can be determined by putting the fusion tree of (i) onto the 
2-D lattice, then smoothly deforming the lattice into a chain 
in accordance with the linearisation prescription, (iii) The 
unitary matrices corresponding to the required F moves and 
braiding operations may be absorbed into M, defining a new 
operator M' on the linearised lattice. 



now assign new tensors (iVCW)t^^ and (7VCCW)t/Ju 
such vertices, such that writing these transformations in 
the notation of Sec. IIV Al is trivial. This is shown in 

Fig.mii). 

Exphcit expressions for the new vertex tensors {N'^^Y 
and (iVCCW)t 

may be obtained by recognising that 
Fig. fTlTi) describes the action of unitary transformations 
on 7V^"^". When the bend is counterclockwise, the cor- 
responding unitary matrix is derived from (A"^)„^, and 
when the bend is clockwise, the unitary matrix is derived 
from {B'^^)uv We will denote these unitary matrices A^^ 
and B^^ respectively. We then have 



to 



7/3 



S ^ 



(13) 
(14) 



and conjugation describes equivalent vertices N'-^^ and 
]SfCCW ^Y^Qj^ a bend is absorbed into a splitting tree. 







(ii) 





i To / 




i To / 



(iii) 





FIG. 11. Vertical bending of legs (i) in the standard dia- 
grammatic notation, and (ii) in the diagrammatic notation of 
Sec. IIV Al White triangles represent Frobenius-Schur indica- 
tor flags, (iii) Legs on the matrix representations of states 
and operators may also absorb bends. 



Knowing how the absorption of bends acts on a vertex 
tensor, we may readily infer how the same process acts on 
the matrix representation of an operator. In Fig. Illf iii) 
we see a bend absorbed into the matrix M^^ , result- 
ing in a new object with two lower multi- indices, M'^j^. 
First we exploit the freedom to introduce fusion with the 
trivial charge (denoted I), with degeneracy 1. The cor- 
responding object takes only one value on its upper 
left multi- index, and is fully defined by N'^lj^^ = 5^^ . Ab- 



sorbing the bend into this fusion vertex as per Eq. p3)) 
yields A^^S^'^S^p — A-^p, which may be then combined 
with M^^ to give 

M'^p = M^-'A^p. (15) 



In conjunction with the relationships given in Fig. I12[ 
this gives us the ability to move a matrix past a bend. 
An example of this is given in Fig. 1131 for which M and 
M' are related according to 



M'^=A^ 



(16) 
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FIG. 12. Opposing pairs of Frobenius-Schur indicators (i) on 
a pair of bends equivalent to the identity, and (ii) on a pair 
of bends such as might be used when computing a quantum 
trace, (iii) As an anyon model can always be specified such 
that the Frobenius-Schur indicators are ±1, reversing a pair of 
contiguous opposed Frobenius-Schur indicator flags is always 
free. 






V 



FIG. 13. Moving a matrix across a bend in a tensor network 
diagram, (i) Initial diagram, (ii) Bends are absorbed into the 
matrix, (iii) New bends are introduced, in accordance with 
Fig. 112^ 1. (iv) A pair of contiguous, opposed Frobenius-Schur 
indicators are reversed, as per Fig. 1 12( 111*). The initial and 
final matrices M and M' are related as specified in Eq. (|16p . 



where represents reversal of the Frobenius-Schur in- 
dicator flag on the lower bend. Finally, bending may 
also allow more efficient contraction of pairs of anyonic 
operators, as shown in Fig. UM 

Having described the action of bends, it is custom- 
ary also to introduce a second type of F move which is 
described by the tensor {F^^^^)(ar,uiU2)iaeU3Ui) (Fig.^. 
This tensor may be derived from (i^a4''"'")(a5«i«2)(a6«3«4) 
by bending, and as with {F^^"-^"-^)(^a,uiU2)(aeU3U4) these 
F moves perform a transformation of the fusion tree, ac- 
companied by the introduction of a unitary matrix which 
can be absorbed into the matrix representation of the 
operator. These unitary matrices correspond to the con- 
secutive application of a bend, an F move of the original 
type, and a second bend whose action is the inverse of 
the first. 



D. Constructing a tensor network 

Now that we have developed a formalism for anyonic 
tensors, we may convert an existing tensor network al- 
gorithm for use with anyons. First, the tensor network 
must be drawn in such a manner that every leg has a 
discernible vertical orientation. Although these orienta- 
tions may be changed during manipulation of the tensor 
network, an initial assignment of upward or downward 
direction is required. Second, all tensors must be repre- 
sented by entirely convex shapes, such as circles or reg- 
ular polygons. For existing tensor network algorithms 
such as MERA and PEPS, this requirement is trivial. 




FIG. 14. The use of bends may permit the more efficient 
contraction of pairs of anyonic operators. In the sequence of 
events marked (i), operator A is first raised to the space of 
three sites then contracted with _B. In sequence (ii) the op- 
erators are instead contracted using bends. For many anyon 
models the latter approach offers a significant computational 
advantage. 



"2 — S (F i^al) (a^^u„u,)(a^,u,.u, j 



FIG. 15. Now that we may bend legs up and down it is 
customary to introduce a further type of F move, derived by 
applying bends to the one presented in Fig. [TJiii). 



However, it is conceivable that future algorithms might 
involve superoperator-type objects whose graphical rep- 
resentations interleave upward- and downward-pointing 
legs. Concavities on these objects may be eliminated 
by replacing some of their upward-pointing legs with 
downward-pointing legs (or vice versa), followed by a 
bend (Fig. I16f i)-(ii)). A similar treatment may be ap- 
plied to any superoperators which arise during manipu- 
lations of the tensor network, introducing a pair of bends 
as in Fig. I12f i) and then absorbing one into the matrix 
representation of the object. 

If working with an anyon model that has non-trivial 
Frobenius-Schur indicators, then indicator flags must be 
applied to all bends. Initial choices are a matter of con- 
venience, and it is frequently possible to assign these in- 
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FIG. 16. Construction of an anyonic tensor corresponding 
to a normal tensor with more than three legs, (i) The origi- 
nal tensor, (ii) If required, any concavities are eliminated by 
introducing bends, (iii) Frobenius-Schur indicators are as- 
signed to the bends, (iv) Directions are assigned to all legs, 
consistent with the rest of the network, (v) Legs are collected 
together into fusion and splitting trees. The central object, 
representing degrees of freedom of the tensor, now has less 
than four legs, (vi) If desired, bends can be re-absorbed into 
the fusion and splitting trees. 



dicators in opposed pairs, as shown in Fig. [T^l If these 
paired indicators are not flipped or are only flipped in 
adjacent opposed pairs during subsequent manipulations 
of the tensor network, then they may frequently be left 
implicit. 

Next, if there exist charges in the anyon model which 
are not self-dual, a direction (represented by a solid ar- 
row) must be assigned to every multi-index. Any tensor 
with more than three legs (e.g. M in Fig. [T6|) is then 
replaced by a trivalent tensor network consisting of a 
core object, e.g. M^^ , which contains the free param- 
eters of the tensor, and as many copies of N or N'^ as 
are required to provide the correct output legs. These 
tensors N, iV^ correspond to vertices in the fusion and 
splitting trees associated with M^'^, yielding the corre- 
sponding anyonic tensor. Objects with three legs or less 
can be directly identified with an anyonic tensor object 
carrying the appropriate number of indices (i.e. three 
multi-indices and a vertex index u), though for consis- 
tency with the methods described in Sections IIVBI and 
llVCI we point out that it is possible to similarly replace 
three-legged objects with anyonic operators consisting of 
a central matrix M^^ and a fusion or splitting vertex, if 
desired. 

Any bends introduced earlier may now be reab- 
sorbed, so that some vertices now correspond to (iV^^), 
(iVCCW)^ (iVCW)t^ and (iVCCW)t. This step, however, 
is optional as it may be more convenient for subsequent 
manipulations of the tensor network if the bends are left 
explicit. The anyonic tensors are then connected pre- 
cisely as in the original ansatz. 

Manipulations of the anyonic tensor network are equiv- 



alent to those performed on the spin version of the ansatz, 
differing only in that the degrees of freedom of the ten- 
sor network are now expressed entirely by the at-most- 
trivalent central objects, and certain topological elements 
such as braids and vertical bends must be accounted for 
in accordance with the prescriptions of Sec. IIV CI These 
changes may naturally imply minor changes to the ma- 
nipulation algorithms, and we will see examples of this 
in the 1-D MERA. Similar considerations will apply to 
other tensor network algorithms. 

Our construction of an anyonic tensor network draws 
upon two important elements which have previously been 
observed in other, simpler, physical systems: 

1. Tensors in the ansatz exhibit a global symmetry, 
which may be nonabelian. Exploiting a nonabelian 
symmetry requires that the ansatz be written in 
the form of a trivalent tensor network. This has 
previously been observed and implemented for non- 
abelian Lie group symmetries such as SU{2)^ 

2. Tensors in the ansatz must be able to account 
for non-trivial exchange statistics. This has pre- 
viously been observed in the simulation of systems 
of fermions^^— where efficient implementation of 
particle statistics can be achieved through the use 
of "swap gates" ^.^^^ 

In both cases, anyonic tensor networks extend the con- 
cepts introduced in previous work. The symmetry struc- 
ture of an anyon model may be a quantum group, for 
example a member of the series SU{2)k, k G Z,+ , rather 
than having to be a Lie group, and this permits repre- 
sentation of nonabelian anyonic systems whose Hilbert 
space does not admit decomposition into a tensor prod- 
uct of local Hilbert spaces. Similarly, anyonic braid- 
ing may be implemented using a generalisation of the 
fermionic "swap gate" formalism. When braiding, parti- 
cle exchange may introduce transformation by a unitary 
matrix rather than by a sign, and efficient implementa- 
tion of the resulting swap gates is particularly important 
for the simulation of 2-D systems. 

Although anyonic systems pose a number of unique 
challenges, we see that these are addressed by develop- 
ments based on existing techniques, and we therefore an- 
ticipate that the resulting generalisations of existing ten- 
sor network ansatze should still be capable of accurately 
representing the states of an anyonic system. 



E. Contraction of anyonic tensor networks 



The techniques described in Sees. HVB] and HVCl {F 
moves, braids, bending of legs, elimination of loops, dia- 
grammatic isotopy, flipping of Frobenius-Schur indicator 
flags, and the use of N'^^^ tensors) suffice to contract any 
network of anyonic tensors written in the form of matrices 
with degeneracy indices, and unlabelled trees. Through 
careful application of these techniques, and avoiding at 
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all times processes which would yield a tensor with more 
than three legs, the matrix representations of any pair of 
contiguous tensors in a network may always be brought 
into conjunction such that their multi-indices can be con- 
tracted in the manner of Eq. (|12p. and any tensor net- 
work may be contracted by means of a sequence of such 
pairwise contractions. 

That a tensor network may represent a system of 
anyons in this way is possible because throughout the 
anyonic tensor network, each value of a degeneracy index 
is associated with a specific labelling of the correspond- 
ing unlabelled tree. Consequently it is always possible 
to fully reconstruct any operation in terms of the more 
verbose representation of Fig. |4l 

An anyonic tensor network is therefore fully specified 
merely by the unlabelled tree (with Frobenius-Schur in- 
dicator flags if required) , and the values and locations of 
the matrix representations of its tensors, written in the 
degeneracy index form. 



V. EXAMPLE: THE 1-D MERA 
A. Construction 

To construct an anyonic MERA for a 1-D lattice with 
n sites, where n satisfies 71 = 2 x 3*^, fc S Z+, we be- 
gin with a "top" tensor on a two-site lattice Ct whose 
matrix representation is of a computationally convenient 
size. (The top tensor is named for its position in the 
usual diagrammatic representation of the MERA, where 
diagrams with open legs at the bottom correspond to 
a ket. For anyons the converse convention applies, and 
consequently in Fig. [TWi) the "top" tensor is ironically 
located at the bottom.) 

To each leg of the top tensor, we now append an isom- 
etry (Fig. [TWii)). The matrix representations of the 
isometrics consist of rectangular blocks, as described in 
Sec. IIIIBl and we choose isometrics whose fusion trees 
have three legs, so as to construct a ternary MERAi^ 
Next, disentanglers are applied above the isometrics. For 
periodic boundary conditions this must be performed in 
a manner which respects the anyonic braiding rules, as 
shown in Fig. [TWiii) . We identify the open legs of the re- 
sulting network as the sites of a lattice Cr^i, and the rows 
of disentanglers and isometrics may be understood as a 
coarse-graining transformation taking a finer-grained lat- 
tice Cr-i into a coarser-grained lattice similar to the 
standard MERA. Note that the geometry of the periodic 
lattice is reflected by the connections of the disentanglers. 
Specifically, whether the outside legs are braided over or 
under the other lattice sites reflects whether the lattice 
closes towards or away from the observer. 

The application of anyonic isometrics and disentan- 
glers is now repeated k times (Fig.fTWiVfiii) corresponds 
to fc = 1), until the ansatz has n legs. The final row of 
isometrics should be chosen such that each of their upper 
legs have the same charges and degeneracies as the sites 




FIG. 17. Construction of a 1-D ternary MERA on a periodic 
lattice from anyonic operators, (i) The "top" tensor, T. (ii) 
Isometries, w. (iii) Disentanglers, u. The fusion tree repre- 
senting an anyonic state (or ket) is usually drawn with the 
lattice sites at the top, so this MERA has been constructed 
"upside down" when compared with the diagrams in Refs. [2^ 
and [27!. This is unimportant, and we could equally well have 
decided to follow the convention usually adopted in tensor 
network algorithms, labelled the tensors in (i)-(iii) by T^, w^, 
and u\ and identified diagram (i)-(iii) as a bra. (iv) Structure 
of a 2-site term in the Hamiltonian, h, or a 2-site reduced den- 
sity matrix, p. (v) Disentanglers and (vi) isometries satisfy 
the relationships uu' = I, ww'' — I. 



of the physical lattice C{), and the open legs above the 
last row of disentanglers are identified with the physical 
lattice. For coarse-grained lattices Ci to the dimen- 
sions of the lattice sites correspond to the lower legs of 
the isometries and are chosen for computational conve- 
nience, subject to the requirement that each charge sector 
is sufficiently large to adequately reproduce the physics 
of the low-energy portion of the Hilbert space. For all 
other legs, their charges and degeneracies are determined 
by requiring consistency with Eq. ([2]). Initial choices of 
which charges to represent on the "top" tensor and on the 
lower legs of the isometries, and with what degeneracies, 
must be guided either by prior knowledge about the phys- 
ical system, or by balancing computational convenience 
against the inclusion of a broad and representative range 
of possible charges. When used in a numerical optimisa- 
tion algorithm, the choice of relative weightings for the 
different charge sectors may often be refined by exami- 
nation of the spectra of the reduced density matrices on 
the coarse-grained lattices, after initial optimisation of 
the tensor network is complete. 

This concludes construction of the MERA for a state 
on a finite, periodic 1-D anyonic lattice. That this ten- 
sor network does represent an anyonic state is easily seen 
by sequentially raising tensors, performing F moves, and 
combining tensors, until the entire network is reduced to 
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a single vector whose length is equal to the dimension 
of the physical Hilbert space, and an associated fusion 
tree. These then represent the state of the system as 
per Eq. Q . The structure of this tensor network closely 
resembles that of the normal MERA, according to the 
identifications given in Fig. 1171 consequently we an- 
ticipate that it will share many of the same properties, 
including the ability to reproduce polynomially decay- 
ing correlators in strongly correlated physical systems. 
Open lattices may also be easily represented by omitting 
the braided disentanglers at the edge of the diagram. 

We also note that in common with the MERA for spins, 
the anyonic MERA may be understood as a quantum 
circuit, although one which carries anyonic charges in its 
wires. Any junction in the fusion/splitting trees may be 
associated with a or tensor, and the entire network 
may be considered as the application of a series of gates to 
a Hilbert space of fixed dimension beginning mostly (or 
entirely, if the top tensor is considered to be the first gate) 
in the vacuum state, with individual gates introducing 
entanglement across some limited number of wires. 



B. Energy minimisation 

The anyonic MERA can be used as a variational ansatz 
to compute the ground state of a local Hamiltonian. 
The Hamiltonian is introduced as a sum over nearest 
neighbour interactions, each term having the form of 
Fig- [TTT iv). and optimisation of the tensor network is 
carried out in the usual manner;^ Also as per usual, 
Hamiltonians involving larger interactions, such as next- 
to-nearest neighbour, can be accommodated by means of 
an initial exact n-into-one coarse-graining of the physical 
lattice. 

As in Ref. [13, optimisation of the MERA then con- 
sists of repeatedly lifting the Hamiltonian from £o to the 
coarse-grained lattices, updating their isometrics and dis- 
entanglers, and lowering the reduced density matrix, or 
the top tensor and its conjugate. When lifting the Hamil- 
tonian or lowerin g th e reduced density matrix, then the 
diagrams in Ref. |27| taken in conjunction with the key 
given in Fig. [17] serve to describe networks of anyonic op- 
erators which, when contracted to a single operator, yield 
the lifted form of the Hamiltonian or lowered form of 
the reduced density matrix respectively. Similarly, when 
optimising disentanglers or isometrics, the diagrams of 
Ref. and the identifications in Fig. [T7] indicate how to 
construct an anyonic operator which constitutes the envi- 
ronment of the anyonic operator being optimised. How- 
ever, once the admissible ranges of charges and degen- 
eracies on each leg have been fixed, the only optimisable 
content of an anyonic operator is its matrix representa- 
tion. Consequently, the fusion and splitting tree contri- 
butions should be evaluated and absorbed into the op- 
erator and its environment, reducing them both to their 
matrix representations, denoted M and E respectively 
(see Fig. [T5)) . If the singular value decomposition of E 




FIG. 18. (i) Anyonic operator E constitutes the environment 
of operator M. Factors arising from the fusion and splitting 
trees should be evaluated and absorbed into matrices E and 
M, following which (ii) matrix E constitutes the environment 
of matrix M. After (iii) updating the matrix M to M' , (iv) 
the fusion and splitting trees of M should be reinstated, the 
numerical factors associated with this process being the in- 
verse of the fusion tree factors previously absorbed into matrix 
M. Frobenius-Schur flags in (i)-(ii) are represented by white 
triangles, and are not to be confused with the black arrows 
which indicate the orientation of lines in the fusion/splitting 
trees. 

is written E — U SW"^ , then the updated matrix con- 
tent M of the anyonic operator being optimised is given 
by —WU\ minimising the value of Tr(i?M) subject to 
the usual constraint for disentanglers and isometrics that 
MMt = I (Fig. [I7i;v)-(vi)). The fusion/splitting tree 
content of the operator can then be restored, along with 
any appropriate numerical factors that may be required. 

As with the standard MERA, the "top" tensor is con- 
structed by diagonalising the total Hamiltonian on the 
most coarse-grained lattice, -fftot on Cr- As Cr is a two- 
site lattice, the total Hamiltonian -fftot is a sum of two 
terms, H12 and H2i. For the translation- invariant any- 
onic MERA, we may formally define H21 in terms of 
H12 as shown in Fig. [111 and the top tensor T (together 
with any factors arising from the chosen normalisation 
scheme) then corresponds to the lowest-energy eigenstate 

of Htot- 



C. Scale invariant MERA 

Having identified the anyonic counterparts of the ten- 
sors of the standard MERA, and described how these 
tensors may be lifted, lowered, and optimised, the algo- 
rithm for the scale-invariant MERA described in Ref. [26| 
may also be implemented for anyonic systems, simply by 
applying the dictionary of Fig. 1 171 and the techniques de- 
scribed in Sec. IV B] As with optimisation of u and w, the 
computation of the top reduced density matrix (which is 
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Hl2 
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H21 



12 12 12 

12 12 




FIG. 19. Definition of H21 in terms of H12, on the most 
coarse-grained lattice (jCt) of the translation invariant peri- 
odic MERA. Lattice jCt is a two-site periodic lattice. 



a descending eigenoperator of the scaling superoperator 
with eigenvalue 1) may be understood as a calculation of 
the matrix component p^^ of the reduced density matrix 
p. The ascending eigenoperators of the scaling superop- 
erator, or local scaling operators of the theory, may also 
be computed in this manner. 



D. Results 

To demonstrate the effectiveness of the anyonic gen- 
eralisation of the MERA, we applied it to a 1-D critical 
system of anyons whose physical properties are already 
well known: The golden chain.** This model consists of a 
string of Fibonacci anyons subject to a local interaction. 
Fibonacci anyons have only two charges, 1 (the vacuum) 
and T, and one non-trivial fusion rule (r x r — > 1 + r). 
The simplest local interactions for a chain of Fibonacci t 
anyons are nearest neighbour interactions favouring fu- 
sion of pairs into either the 1 channel (termed antifer- 
romagnetic, or AFM), or the r channel (termed ferro- 
magnetic, or FM). Both choices correspond to critical 
Hamiltonians, associated with the conformal field the- 
ories A^(4, 3) and M{5,4) for AFM and FM couplings 
respectively. Individual lattice sites are each associated 
with a charge of r. 

The AFM and FM Hamiltonians act on pairs of ad- 
jacent Fibonacci anyons. On a pair of lattice sites each 
carrying a charge of r, the matrix representations of the 
AFM and FM Hamiltonians are written 




FIG. 20. Determination of eigenoperators {(p) and associated 
scaling dimensions (A) for the one-site scaling superoperator 
of the anyonic 1-D MERA. Eigenoperators may be classified 
according to the charges on edges yi and 2/2. One interpreta- 
tion of these labels is that, in addition to the sites of the 1-D 
lattice, there may exist free charges lying in front of and be- 
hind the anyon chain. The labels j/i and 1/2 then represent the 
transfer of charge between these regions and the 1-D lattice. 



where a multi-index value of 1 corresponds to the vac- 
uum charge, 2 corresponds to r, and the charges are 
non-degenerate. We optimised a scale-invariant MERA 
on the golden chain for each of these Hamiltonians, and 
computed local scaling operators using the tensor net- 
work given in Fig. [201 The operators calculated using 
this diagram may be classified according to the values 
of the charge labels yi and y2, and the scaling dimen- 
sions and conformal spins which we obtained are given 
in Tables m] and HVl and Fig.mi 

Comparison of the AFM case with existing results in 
the literature show that the scaling dimensions obtained 
when yi = y2 correspond to those obtained when study- 
ing a system of anyons with a toroidal fusion diagrami^ 
For a system of anyons on the torus it is possible to de- 
fine an additional topological symmetry* and classify lo- 
cal scaling operators according to whether or not they 
respect this symmetry. Operators satisfying yi — y2 — ^ 
correspond to those which respect the topological sym- 
metry, and those satisfying yi = y2 = t do not. We 
will discuss the interpretation of the different sectors and 
their relationship to anyons on the torus in a forthcoming 
paper.— 



When yi y2 the scaling operators obtained are chi- 
ral, with those obtained from yi — 1, y2 = t and 
yi — T, 2/2 = 1 believed to form conjugate pairs. 
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yi = 2/2 = 


1 




yi = y2 = r 


Exact 


Numerics 


Error 


Exact 


Numerics 


Error 








0% 




3/40 


0.0751 


+0.19% 


7/8 


0.8995 


+2.80% 




1/5 


0.2006 


+0.28% 


7/8 + 1 


1.9096 


+1.85% 




3/40 + 1 


1.0730 


-0.19% 


7/8 + 1 


1.9141 


+2.09% 




3/40 + 1 


1.0884 


+ 1.25% 


+ 2 


2.0124 


+0.62% 




6/5 


1.2026 


+0.21% 


+ 2 


2.0181 


+0.90% 




1/5 + 1 


1.2156 


+1.30% 



(i) 





= 1, y2 = 


T 




yi 


= r,y2 = 


1 


Exact 


Numerics 


Error 


Exact 


Numerics 


Error 


19/40 


0.4757 


+0.14% 




19/40 


0.4757 


+0.14% 


3/5 


0.6009 


+0.15% 




3/5 


0.6009 


+0.15% 


19/40+1 


1.4549 


-1.37% 




19/40+1 


1.4549 


-1.37% 


19/40+1 


1.5022 


+1.85% 




19/40+1 


1.5022 


+1.85% 


3/5 + 1 


1.5414 


-3.66% 




3/5 + 1 


1.5414 


-3.66% 


3/5 + 1 


1.6129 


+0.80% 




3/5 + 1 


1.6129 


+0.80% 



TABLE III. Scaling dimensions for Fibonacci anyons with an- 
tiferromagnetic nearest neigbour interactions on an infinite 
chain. Numerical values were computed using an anyonic 
MERA with maximum degeneracies for charges 1 and r of 
3 and 5 respectively (denoted x ~ [3,5]), and are grouped 
according to their classification by the values of y\ and j/2 in 
Fig. [201 



yi = 2/2 = 1 




2/1 = y2 = r 


Exact 


Numerics 


Error 




Exact 


Numerics 


Error 








0% 




2/15 


0.1329 


-0.35% 


4/3 


1.3514 


+1.36% 




2/15 


0.1339 


+0.44% 


4/3 


1.3695 


+2.71% 




4/5 


0.8134 


+ 1.67% 


+ 2 


1.9519 


-2.41% 




2/15 + 1 


1.0937 


-3.49% 


+ 2 


1.9742 


-1.29% 




2/15 + 1 


1.1108 


-1.99% 


1 + 4/3 


2.2570 


-3.27% 




2/15 + 1 


1.1622 


+2.55% 
















= 1, J/2 = 


T 




2/1 


= T, 2/2 = 


1 


Exact 


Numerics 


Error 




Exact 


Numerics 


Error 


2/5 


0.3993 


-0.18% 




2/5 


0.3993 


-0.18% 


11/15 


0.7327 


-0.09% 




11/15 


0.7327 


-0.09% 


11/15 


0.7392 


+0.80% 




11/15 


0.7392 


+0.80% 


2/5 + 1 


1.3699 


-2.15% 




2/5 + 1 


1.3699 


-2.15% 


2/5 + 1 


1.3823 


-1.26% 




2/5 + 1 


1.3823 


-1.26% 


11/15+1 


1.6450 


-5.10% 




11/15+1 


1.6450 


-5.10% 



TABLE IV. Scaling dimensions for Fibonacci anyons with fer- 
romagnetic nearest neigbour interactions on an infinite chain. 
Numerical values were computed using an anyonic MERA 
with maximum degeneracies for charges 1 and r of 3 and 5 
respectively (denoted x = [3,5]), and are grouped according 
to their classification by the values of y\ and 1/2 in Fig. 1201 
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FIG. 21. COLOR ONLINE. Scaling dimensions of leading 
primary operators and their descendants, computed for (i) 
antiferromagnetic and (ii) ferromagnetic local Hamiltonians 
(|17|) on the golden chain. Results are grouped into conformal 
towers, with a slight horizontal spread introduced to show 
the degeneracies of the descendant fields. A circled cross indi- 
cates a primary field, and a plain cross indicates a descendant. 
Dashed lines indicate values predicted from CFT. 



VI. SUMMARY 

Numerical study of systems of interacting anyons is 
difficult due to their non-trivial exchange statistics. To 
date, study of these systems has been restricted to ex- 
act diagonalisation, Matrix Product States (MPS) for 1- 
D systems, or special-case mappings to equivalent spin 
chains. This paper shovifs hovi^ any tensor netvifork ansatz 
may be translated into a form applicable to systems of 
anyons, opening the door for the study of large sys- 
tems of interacting anyons in both one and two dimen- 
sions. As an example, this paper demonstrates how the 
MERA may be implemented for a 1-D anyonic system. 
This ansatz is particularly important as many 1-D sys- 
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terns of anyons are known which exhibit extended critical 
phases. The structure of the MERA is known to be 
particularly well suited to reproducing long range corre- 
lations, and the scale-invariant MERA has the additional 
advantage of providing simple and direct means of com- 
puting the scaling dimensions and matrix representations 
of local scaling operators. 

We applied the scale invariant MERA to infinite chains 
of Fibonacci anyons under antiferromagnetic and ferro- 
magnetic nearest neighbour couplings, and identified a 
large number of local scaling operators. Our results for 
the scaling dimensions are in agreement with those pre- 
viously obtained by exact diagonalisation of closely re- 
lated systems, and for the relevant primary fields they 
are within 2.8% of the theoretical values obtained from 
conformal field theory. We thus demonstrate that an 
anyonic MERA with x — [3, 5] permits conclusive identi- 
fication of the relevant conformal field theory, and gives a 
level of accuracy comparable to that of the scale invariant 



MERA on a spin chainj2^ 

The anyonic generalisation of the 1-D MERA pre- 
sented here is useful in its own right, but the great- 
est significance of the approach described is that it is 
equally applicable to 2-D tensor network ansatze, and 
hence opens the door to studying the collective behaviour 
of large systems of anyons in two dimensions by numer- 
ical means, in situations where analytical solutions may 
not be possible. 

The authors acknowledge the support of the Australian 
Research Council (FF0668731, DP0878830, DP1092513, 
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Perimeter Institute for Theoretical Physics. 

Note — While preparing this pape r for publication, w e 
became aware of related work by ^ Konig and Bilginl 4^ 
They also present the anyonic 1-D MERA, providing 
proof of principle by computing ground state energies for 
finite systems of Fibonacci anyons with x— [Ij 1] [s — 2 
in their notation). 



* |pfeifer@physics. uq.edu.au] 

^ A. Y. Kitaev, Annals of Physics 303, 2 (2003). 

^ C. Nayak, S. H. Simon, A. Stern, M. Freedman, and 



S. Das Sarma, plev Mod. Phys. 80, 1083 (2008) 
J. S. Xia, W. Pan, C. L. Vicente, E. D. Adams, N. S. Sulli- 
van, H. L. Stormer, D. C. Tsui, L. N. PfeifTer, K. W. Bald- 
win, and K. W. West, Phys. Rev. Lett. 93, 176809 (2004) 
A. Feiguin, S. Trebst, A. W. W. Ludwig, M. Troyer, 
A. Kitaev, Z. Wang, and M. H. Freedman, 
iPhys. Rev. Lett. 98, 160409 (2007), 



S. R. White, |Phys. Rev. L ett. 69, 2863 (1 992)| 

G. Sierra and T. Nishino, Nuc. Phys. B 495, 505 (1997)] 
W. Tatsuaki, Phys. Rev. E 61, 3199 (2000) 
S. Trebst, E. Ardonne, A. Feiguin, D. A. 
Huse, A. W. W. Ludwig, and M. Troyer, 



Phys. Rev. Lett. 101, 050401 (2008)| 



F. C. Alcaraz, M. N. Barber, M. T. Batchelor, R. J. Baxter, 
and G. R. W. Quispe l, J. Phys. A 20, 6397 (1987). 
S. Todo and K. Kato, |Phys. Rev. Le tt. 87, 047203 (2001) 
F. Verstraete and J. I. Cirac, 

ajXivxond^-mat/0407066vl (2004)[ 

T. Nishino and K. Okunishi, 

J. Phys. Soc. Jap 67, 3066 (1998)| 

and 



Z.-C. 



Gu, 



M. 



Levin, 



X.-G. 



Wen, 



Phys. Rev. B 78 , 205116 (2008) 

Z. Y. Weng, 



H 



Jiang, 



and T. Xiang, 



Phys. Rev. Lett. 101, 090603 (2008) 



J. Jor dan, R. Orus, G. Vidal, F. V erstraete, 
Cirac, 'Phys. Rev . Lett. 10 1, 25060 2 (2008)| 



and J. I. 



Tagliacozzo, G. Evenbly, 
Phys Re v. B 80, 235127 (2009) 
L! Cincio, J. Dziarmaga, 



and G. Vidal, 



and M. M. Rams, 



Phys. Rev. Lett. 100, 240603 (2008) 
H 



Evenbly 



and 



Phys. Rev. Lett. 102, 180406 (2009) [ 



G. Evenbly and 

Phys. Rev. Lett. 104, 187203 (2010)] 



Vidal, 



Vidal, 



' P. Di Francesco, P. Mathieu, and D. Senechal, Conformal 

field theory (Springer, New York, 1997). 

G. Vidal, Phys. Re v. Lett. 99, 220405 (2007) 

G. Vidal, Phys. R ev. Lett. 101, 110501 (2008) 
' G. Vidal, in Understanding Quantum Phase Transitions, 

edited by L. D. Carr (Taylor & Francis, Boca Raton, 2010). 

X. Chen, Z.-C. Gu, and X.-G. Wen, 

ar Xiv:1004.3835vl [con d-mat.str -el] (201 0) 

V. Giovannetti, S. Montang ero, and R. Fazio, 

|Phys. Rev. Lett. 101, 18050 3 (2008)', 

' R. N. C. Pfeifer, G. Evenbly, and G. Vidal, 
'Phys. Rev. A 79, 040301 (2009) 



G. Evenbly and G. Vidal, Phys. Rev. B 79, 14 4108 (2009) 
A. Kitaev, Annals of Physics 321, 2 (2006) 
P. H. Bonderson, Non-Abelian Anyons and Interferometry, 
Ph.D. thesis, California Institute of Technology (2007). 
S. Singh, R. N. C. Pfeifer, an d G. Vidal, 
|arXiv:0907.2994vl [cond-mat.str-el] (2009) 



P. Corboz, G. Evenbly, F. Verstraete, and G. Vidal, 
Phys. Rev. A 81, 010303 (2010) 

C. V. K raus, N. S chuch, F. Verstraete, and J. I. Cirac, 
Phys. Rev. A 81, 052338 (2010), 

and 



Pineda, 



Barthel, 



Eisert, 



Phys. Rev. A 81, 050303 (2010) 

P. Corboz, R. Orus, B. Bauer, and G. Vidal, Phys. Rev. 
B 81, 165104 (2010). 

T. Barthel, C. Pineda, and 



Eisert, 



Phys. Rev. A 80, 042333 (2009) | 



Q.-Q. Shi, S.-H. Li, J.-H. Zhao, and H.-Q. 
arXiv:0907.5520vl [cond-ma t.str-e l] (2009)| 



Zhou, 



I. Pizorn and 

Phys. Rev. B 81, 245110 (2010) 
ZA:. Gu,"" F. Verstraete, 



Verstraete, 



and X.-G. Wen, 



: arXiv:1004.2563vl [cond-mat.str-el] (2010), 
R. N. C. Pfeifer et a/., |in preparation! 
S. Trebst, M. Troyer, Z. Wang, and A. W. W. Ludwig, 
Prog. Theo. Phys. Supp. 176, 384 (2008), 



R. Konig and E. Bilgin, ,Phys. Rev. B 82, 125118 (2010) 



